Structural and electronic properties of Si/Ge nanoparticles 

Abu Md. Asaduzzaman* and Michael Springborg^ 
Physical and Theoretical Chemistry, 
University of Saarland, D-66123 Saarbriicken, Germany 
(Dated: February 6, 2008) 

Abstract 

Results of a theoretical study of the electronic properties of (Si)Ge and (Ge)Si core/shell nanopar- 
ticles, homogeneous SiGe clusters, and Ge|Si clusters with an interphase separating the Si and Ge 
atoms are presented. In general, (Si)Ge particles are more stable than (Ge)Si ones, and SiGe sys- 
tems are more stable than Ge|Si ones. It is found that the frontier orbitals, that dictate the optical 
properties, are localized to the surface, meaning that saturating dangling bonds on the surface 
with ligands may influence the optical properties significantly. In the central parts we identify a 
weak tendency for the Si atoms to accept electrons, whereas Ge atoms donate electrons. 

PACS numbers: 73.22.-f, 61.46.Df, 73.61.Le, 36.40.-c 



* e-mail: a.asaduzzaman@mx.uni-saarland.de 
t e-mail: m.springborg@mx.uni-saarland.de 



1 



I. INTRODUCTION 



Traditionally, materials properties have been controlled by varying structure and com- 
position of the materials. During the last quarter of a century a new parameter has been 
added, i.e., size. The fact that the materials properties change drastically when the dimen- 
sion of the materials becomes comparable with the typical length scale of the phenomenon 
of interest together with the ability to control the production of materials in this size range 
has led to the development of 'nano-science'. 

One of the materials classes where materials in the nm range are expected to have a large 
impact on the development of new and/or better devices is semiconductors. For those, mod- 
ified electronic properties may show up when the size of the nanocrystals is comparable with 
the spatial extension of the excitons. For instance, the light emission properties of semicon- 
ductor nanocrystal quantum dots, specifically the tuning of color afforded by the quantum 
size effect, is important for application of these materials in light emitting device o 1 ' 2 ! 3 ' 4 and 
as biological fluorescence markers j^ 4 ^ 

A further development which has opened up new possibilities to control and vary the 
materials properties is the successful production of core/shell nanoparticles, consisting of a 
core of one material coated by a well-defined shell layer of another material. These systems 
exhibit unique and advanced properties over single- component nanoparticles, making them 
attractive for use in a wide range of real-world applications, and are therefore of extensive 
scientific and technological interest. In the last few years, much effort has been focused on the 
synthesis, fabrication, and characterization of the core/shell structured semiconductor het- 
erostructures with tailored properties. The growth of the shell on the core material to form a 
core/shell heterostructures has been successfully demonstrated on the surface reconstruction 
of nanostructured material. Among the ingredients that dictate the electronic, electrical, 
optical, and chemical properties of core/shell nanostructures, the surface-to- volume ratio, 
the shell type and shell thickness are important, although a precise understanding of the 
relations between structure, size, and composition on the one hand and property on the 
other hand is lacking. 

Recently, Lauhon et ai.— reported the epitaxial growth of crystalline silicon-germanium 
and germanium-silicon core/sheath structures. Despite a 4% lattice mismatch of the macro- 
scopic crystalline systems, the experimental study demonstrated that for Si core nanowire, 
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the Ge shell is fully crystallized at low temperatures and for a Ge core nanowire, an amor- 
phous Si shell is formed initially and after thermal annealing the shell becomes crystallized. 
In another study, Malachias et al.— experimentally showed the formation of Ge domes with a 
Si core and a Ge shell. Kolobov et al— in their experimental study showed the formation of 
nanocrystals of a Ge core with a SiGe shell. On the other hand, there are few experimental 
and theoretical studies on core/shell materials reported . 11 ^ 2 ' 13 Most studies on core/shell 
studies have considered binary compounds like CdSe/ZnS, ZnS/CdS, CdSe/CdS, etc. Core- 
shell studies on pure elements are very scarce. Recently, Musin and Wangi^ reported a 
theoretical study on the epitaxial Si-Ge core/shell structure, a theoretical realization of the 
experimental work of Lauhon et al. 

In this paper, we present the results of a theoretical study of structural and electronic 
properties of naked Si-Ge and Ge-Si core/shell nanoparticles. Musin and Wang took into 
consideration the epitaxial growth of Si-Ge core/sheath nanowires and studied the compo- 
sitional dependency of the structural parameters and the band gap energy. Here, we have 
considered naked Si-Ge and Ge-Si core/shell nanoparticles for which we assumed that the 
structure is related to that of a spherical cut-out of the infinite crystal with a zincblende- 
or diamond-like structure. Although our assumption may affect the results of the calcula- 
tions, we believe that by considering a larger number of sizes (almost 100 sizes with up to in 
total almost 200 atoms) our study allows for drawing general conclusions on those systems. 
Starting with a small core of only 8 atoms and a thin shell of 24 atoms we have gradually in- 
creased the sizes of both core and shell. In order to find out how the properties of core/shell 
particles depend on the two elements, we have, in addition, carried out calculations on pure 
Si clusters, pure Ge clusters, homogeneous SiGe clusters, and spherical structures with half 
of the sphere made up of Si and the other half made up of Ge atoms (Si|Ge). A representa- 
tive example of a core/shell particle as well as of homogeneous SiGe and of a Si|Ge particle 
is depicted in Fig. ^ We mention that the experimentally studied systems are considerably 
larger (containing up to several 1000s of atoms), which is very difficult to study theoretically. 
Therefore, in this study we have limited ourselves to a detailed study of different core/shell 
structures along with pure Si and Ge, homogeneous SiGe, and Si|Ge systems with up to 
around 200 atoms, which is around 2 nm in radii. We have optimized all the structures to 
their nearest local total-energy minima whereby all atoms were allowed to move. 
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II. COMPUTATIONAL OUTLINE 



We have used a parameterized density-functional tight-binding method that has been 
described in detail elsewhere . 15 i 16 i 17 The total energy relative to the isolated atoms of a 
given system is written by 

E b = £ * - E ^ + \ £ - £l) (!) 

Here, is the energy of the zth orbital for the system of interest and 6jf- is the energy of the 
jth orbital for the isolated kth atom. Um is a pair potential between the kth and Zth atom. 
The valence single-particle eigenfunctions ipi{r) to the Kohn-Sham equation are expanded 
in a set of atom-centered basis functions Xkim{^), where k denotes the atom and (l,m) the 
angular dependence. The effective one-electron potential in the Kohn-Sham Hamiltonian 
is approximated as a superposition of the atomic potentials of the corresponding neutral 
atoms, and it is assumed that the matrix elements (Xk 1 i 1 m 1 \Vj\xk 2 i2m 2 ) (with Vj being the 
atomic potential at atom j) are vanishing unless k\ = j or k2 — j. Thus, only two-center 
Hamiltonian matrix elements are considered and calculated exactly within the Kohn-Sham 
basis for the diatomic molecules. Finally, the pair potentials U^i are determined so that the 
binding energy curve of the diatomics are well reproduced. Only the 3s and 3p functions of 
Si and the 4s and Ap functions of Ge were explicitly included in the calculations, whereas 
all other electrons were treated within a frozen-core approximation. 

It is obvious that, the approach we are using has been designed for the smallest possible 
systems Si2, Ge2, and SiGe. We have verified the capability of this method to do calculations 
for larger systems by studying some characteristics of bulk Si and Ge. The experimental 
lattice constants of crystalline Si and Ge are 5.43 and 5.66 A, respectively, whereas our 
calculations give 5.46 and 5.71 A, respectively, i.e., within less than 1% of the experimental 
values. The experimental band gaps of bulk Si and Ge are 1.12 and 0.66 eV, respectively, 
and our calculated values are 1.097 and 0.65 eV, respectively which are also close to the 
experimental values. Here, the standard problem of density-functional calculations to yield 
too small band gaps seems to be absent, mainly due to the fact that our basis set is minimal 
in size. 

In addition to the clusters containing both Si and Ge atoms, we have also studied small 
clusters of pure Si or Ge, separately. Among the different arrangements of three atomic Si 
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and Ge clusters, clusters with G% v symmetry are those of the minimum energy configuration. 
Clusters with Dih symmetry are the minimum energy structures among the four atomic 
clusters for both Si and Ge. Earlier theoretical studies on Si and Ge^ and on Si^S also 
reported those minimum energy structures. 

In all cases we relaxed initial structures that we constructed by cutting out a spherical 
part of a diamond-like crystal. The initial lattice constant was taken as that of the pure 
crystalline systems (for the pure Si and Ge clusters) or as the average of those two values 
(for the heteroatomic structures). The center of the sphere was taken as the midpoint of a 
nearest-neighbor bond, giving that the number of atoms would be 2, 8, 20, 32, 38, 56, 74, 
86, 116, 130, 166, or 190 if 1, 2, . . ., 12 atomic shells were included in the initial structure 
(here, we define an atomic shell as being the set of atoms that has the same distance to the 
center of the spherical cut-out in the initial structure). Subsequently, the initial structure 
was allowed to relax to its closest total-energy minimum, whereby all atoms (i.e., both in 
the inner part and in the surface region) were displaced until the forces on them vanish. 



III. RESULTS 

We studied in total 95 different structures. Each of those structures consists of iVg£ $ Si 
atoms and N Gc ^ Ge atoms, i = 1, 2, . . . , 95. Using a least-squares fit we approximated the 
binding energy of those 95 structures by a sum of atomic energies, 

E bfl ~ E Si N Si>i + E Ge N GCii = E b>i . (2) 

Subsequently, we defined one stability energy for each cluster, 

AEi = E b;i — E b;i (3) 

which is the more negative the more stable the cluster is. Finally, we analyse this quantity 
per atom, i.e., 

AEx/N = AE 1 /(iV SM + N Geii ). (4) 
We also considered the stability energy 

AE 2 /N = E b4 /(N Si4 + N Ge!i ). (5) 

Each of those quantities is analysed as a function of the number atomic shells either in 
the core, in the shell, or in the complete nanostructure. I.e., with N t being the total number 
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of atomic shells, 1 < Nt < 12, for the core/shell particles we have N c atomic shells in the 
core and N s = N t — N c shells in the shell part , whereas N s = N t , N c = for the homogeneous 
SiGe clusters and for the Si|Ge nanoparticles. 

The fit of Eq. © resulted in E Si = -2.37 eV and E Ge = -3.58 eV. The fact that E Ge is 
more negative than E$i implies that it is energetically more favorable for Ge atoms than for 
Si atoms to be incorporated into those nanostructures, although the difference in the two 
energies is relatively small. On the other hand, the cohesive energy of the solids equals 4.63 
and 3.85 eV/atom for Si and Ge, respectively (see, e.g., giving that for the elemental 

solids it is energetically more favorable for Si than for Ge atoms to be incorporated into the 
solids. The difference between our results and those for the elemental solids may be due to 
the differences in the systems (both with respect to composition and regarding size) and in 
the finite number of structures that we include in our fit. 

The fact that E Ge < E^i also means that when comparing AEi and AE 2 , the larger the 
number of Ge atoms is compared with that of Si atoms, the lower is AE 2 compared with 
AE 1 . 

In Figs. El and El we show the two energies AE\/N and AE 2 /N as a function of the number 
of shells in either the core region, in the shell region, or in the complete cluster. 

First, the results for the Si|Ge and the homogeneous SiGe nanosystems show clearly that 
the latter is more stable than the former. In addition, in both cases the binding energy 
per atom is a decaying function of the size of the system. This behavior is often found for 
clusters, and can be related to the presence of the surface atoms for which the appearance 
of dangling bonds and lower coordinations lead to energetically less favorable situations. 
When the system size is increased, the relative importance of the surface is reduced, leading 
to an overall decaying total energy per atom as a function of size. The fact that the precise 
number of surface atoms is not completely regular as a function of cluster size can explain 
the deviations from a completely smooth behavior. 

Comparing Figs. Efb) andEfb) the finding that the latter is fairly unstructured, whereas 
the former shows a clear decaing behavior as function of size suggests that for the (Ge)Si 
systems, i.e., a Ge core covered with Si shells, the total energy for this system depends only 
marginally on its size. 

On the other hand, Figs. Efd) shows that for the (Ge)Si nanoparticles the variation in 
the total energy is mainly determined by the number of atomic Si shells in the shell part 
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and largely independent of the number of atomic Ge shells in the core region. A similar 
result is found for the (Si)Ge nanoparticles in Fig. Etc) with, however, some few exceptions 
that are recovered in all results for these systems, i.e., in Figs. I2fa,c,e) and Figs. E(a,c,e). 
It turns out that these exceptions are among the smallest systems we have studied with 
only two atomic shells of Si covered with a low number of Ge shells. We mention that the 
systems with no atomic Si shells of Si, are less stable, so that the (Si)Ge nanoparticles with 
two atomic shells of Si covered with a small number of Ge shells indeed are the most stable 
systems of the present study. We do not have a precise explanation for this finding. 

A further result can be seen when comparing the (Si)Ge core/shell systems with the 
(Ge)Si ones: the former are more stable than the latter. Since the (Si)Ge systems contain 
a Ge surface, whereas the (Ge)Si systems contain a surface of Si atoms, this result may be 
a simple consequence of the lower surface energy of Ge compared with that of Si (see, e.g., 



2l|). 



The HOMO-LUMO gap, E gap , (i.e., the energy gap between the highest occupied molec- 
ular orbital and the lowest unoccupied molecular orbital) is shown in Fig. 0] similarly to 
Figs. El and El With no exception the gap is significantly smaller than those of the elemental 
solids and in many cases it is even close to vanishing. The systems with the largest gaps 
are the core/shell nanoparticles with the smallest number of atomic shells in the complete 
system as well as the Si|Ge systems. Marsen et ai. , 23 ' 24 in their experimental study, showed 
that the band gap of small silicon clusters is larger (the largest value being 0.45 eV) than 
that of the large clusters (most clusters have gap ranging from 0.025 to 0.10 eV) and the 
gap is size dependent. Besides that, they observed essentially no correlation between size of 
the gap and size of the cluster, which is in good agreement with our results. 

In an earlier work on a number of stoichiometric II- VI and III-V nanoparticles we found 
a close correlation between the HOMO-LUMO gap and the stability. 22 In order to explore 
whether a similar correlation exists here, we show in Fig. El-Egap as a function of AEx/N. It 
is clear that no correlation is found in the present case. Replacing AEi/N by AE2/N does 
not change this conclusion. 

Ge/Si nanowires with a Si sheath covering a Ge core have been of some interest recently 
(see, e.g., a l so as active components in semiconductor devices. In these, it is assumed 

that a hole gas is formed in the Ge wire. If the systems were isolated, this would imply a net 
electron transfer from Ge to Si, which, then, also may occur for the system of the present 
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study. In order to study this, we first determine the center of the cluster with n Si atoms 
and m Ge atoms, 

i n+m ^ 

R = — — J2 Rv (6) 

n + m p{ 

and, subsequently, for each atom its so-called radial distance 

r 3 = \R 3 -R \. (7) 

Subsequently, we plot the Mulliken gross populations of the individual atoms as a function 
of the radial distance; cf. Fig. H3 Indeed, it is seen that for the (Si)Ge core/shell particles 
there is a slight tendency for the atomic populations on the Si atoms to be larger than 4, 
whereas for the (Ge)Si core/shell particles, the Ge atomic populations are on the average 
below 4. On the other hand, in the shell region, the trend is less clear. However, in total, 
these core/shell systems do show tendency towards electron- or hole- gas formation in the 
central part, depending on whether Si or Ge is the material of the core. The finding that 
the atomic populations deviate most from the value of 4 for the atoms closest to the surface 
is also observed for the Si|Ge systems. For the homogeneous SiGe systems we also observe 
a small electron transfer from Ge to Si in the central part, which, however, not is the case 
for the central parts of the Ge|Si systems. 

Much of the interest in semiconductor nanoparticles is connected with their partly con- 
trollable optical properties that first of all are dictated by the properties of excitons, which 
in turn are determined by the energy and the charge distribution of the orbitals closest to the 
Fermi level, i.e., the HOMO and the LUMO. Optical relaxation processes strongly depend 
on the spatial distribution of these orbitals. Therefore, we shall study these. To this end we 
construct for any orbital the density 

/2n>\ 3 / 2 

M^l ex P [-a(f-R 3 n (8) 

Here, N i3 is the Mulliken gross population for the jth atom and ith orbital, and a is chosen 
'reasonably', so that illustrative figures result. 

Fig. shows the resulting schematic representation of this radial dependence of the 
HOMO and LUMO. It is seen that in most cases, both the HOMO and the LUMO are lo- 
calized to the surface of the clusters, irrespectively of whether we consider (Si)Ge or (Ge)Si 
systems. In an earlier study on CdSe/CdS core/shell nanoparticles^ it was found that one 
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may find systems for which the LUMO and the HOMO were localized in different parts (i.e., 
in the shell and in the core) of the system, but this is obviously not the case for the present 
systems. Since the gap of crystalline Ge is smaller than that of crystalline Si, and since 
the HOMO and LUMO of the present systems are located to the surface region, it may be 
suggested that the gap of the (Si)Ge particles is smaller than that of the (Ge)Si systems. 
Actually, such a difference has been found for the pure Ge and Si systems by Melnikov and 
Chelikowsky. 27 However, in the present study such a trend is only very marginally found 
for the core/shell systems, which is in agreement with the theoretical results of Musin and 
WangiA who studied Ge/Si and Si/Ge core/sheath nanowires. Moreover, Musin and Wang 
found a larger band gap than those of the pure, crystalline elements for all systems, which 
most likely is due to the passivation of the surface included in their study. This once again 
indicates that without passivation the frontier orbitals are localized to the surface region. 

Also for the homogeneous SiGe systems the frontier orbitals are located to the surface, 
as can be seen in Fig.lHl For the Si|Ge systems we find also this behavior, cf. Fig. El In both 
cases we find that the HOMO is localized to both Si and Ge atoms, whereas the LUMO is 
localized mainly to the Ge atoms. 

IV. CONCLUSION 

In this paper, we have presented the results of our study of the structural and electronic 
properties of naked (Si)Ge and (Ge)Si core/shell nanoparticles together with those of pure 
Si and Ge, homogeneous SiGe and Si|Ge systems with a diameter of up to around 2 nm. 
The interesting properties critically depend on the size of core and shell and also on the type 
of core and shell atoms. Although we did not enter a detailed discussion of this issue, we 
emphasize that our results on pure Si and Ge clusters, in particular concerning the reduced 
band gap compared with the crystalline material, are in excellent agreement with results of 
other experimental and theoretical studies. 

Due to difference in surface energy, (Si)Ge core/shell systems are more stable than are 
(Ge)Si systems. Moreover, the stability of the (Si)Ge systems was quite well described in 
terms of the quantity AE\ as a function of the number of core shells, whereas AE2 as a 
function of the number of core shells gave a good description of the stability of the (Ge)Si 
systems. 
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We did observe a marginal tendency of a Ge— >Si electron transfer, irrespectively of which 
system was forming the core and which the shell. A similar effect was not found for the 
homogeneous SiGe systems, but for the interface Si|Ge systems. 

In all systems, the (unpassivated) surface was dictating the properties of the frontier or- 
bitals. The latter were localized to the surface, and their energies occurred in the energy gap 
of the pure, infinite crystals, leading to quite low band gaps. Therefore, surface passivation 
could be a useful means of tuning the optical properties of these systems. In contrast to our 
earlier results on naked AB semiconductor clusters, we do for the present ones not observe 
a correlation between stability and band gap. Finally, for none of the systems we found a 
charge separation upon electronic excitation. 
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FIG. 1: Ball and stick representation of a representative core/shell structure (upper panel). More 
and less dark parts represent the core and the shell, respectively. Both the sideview (left panel) 
and the cross-sectional view (right panel) are presented. Si|Ge (left panel) and homogeneous SiGe 
(right panel) systems are represented in the lower panel. 
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FIG. 2: The variation in the stability energy per atom AE±/N as a function of the number of 
atomic shells for the different systems, as indicated in the panel. iV s heii is the total number of 
shells, N t , in (a), (b), (g), and (h), the number of shells in the core, N c , in (e) and (f), and the 
number of shells in the shell N s , in (c) and (d). Finally, the lines in (c)-(f) connect the values for 
the systems with the same number of atomic shells in (c,d) the core or in (e,f) the shell. 
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FIG. 4: As Fig. but for the HOMO-LUMO energy gap, E i 
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FIG. 5: The H0M0-LUM0 energy gap as a function of the stability energy per atom AE\/N for 
(a,b) the core/shell systems with (a) Ge covering Si and (b) Si covering Ge as well as for (d) the 
homogeneous GeSi systems and (c) the Ge|Si systems. 
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FIG. 7: Schematic presentation of the radial electron distribution of the HOMO (curves pointing 

upward) and the LUMO (curves pointing downward) for the same systems as in Fig. |HI 
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FIG. 8: Schematic presentation of the radial electron distribution of the HOMO (curves pointing 
upward) and the LUMO (curves pointing downward) for homogeneous SiGe nanoparticles with (a) 
32, (b) 56, (c) 74, (d) 86, (e) 116, (f) 130, (g) 166, and (h) 190 atoms. 
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FIG. 9: Schematic presentation of the radial electron distribution of the HOMO (curves pointing 
upward) and the LUMO (curves pointing downward) for the same number of atoms as Fig. |H] but 
for Si|Ge nanoparticles. 
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